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The data on the dispersion of the magnetic excitations of 
CuGeC>3 in the spin-Peierls D-phase are analyzed. On the 
basis of the lattice structure it is shown that even along the 
chains the d = 2 character cannot be neglected. The symme- 
try of the dispersion differs from the one assumed so far. The 
magnetic resonance data is reinterpreted. The possibility of 
interchain rather than intrachain frustration is discussed. 
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The substance CuGeC>3 is the first inorganic example 
of a spin-Peierls system Q . A second compound which is 
intensively investigated is a'-NaV205 |^,|| An inorganic 
sample is particularly interesting since it allows to employ 
all the experimental tools (for a review, see Q). Thus 
the investigation of CuGeC>3 has become a very active 
field both in theory and in experiment. 

A spin-Peierls system is an antiferromagnetic spin 
system coupled to phonons such that below a certain 
temperature Tsp the translational symmetry is sponta- 
neously broken. The system becomes dimerized due to 
a static lattice distortion which renders the spin cou- 
pling alternating along a certain axis. There are stronger 
bonds (J) and weaker bonds (A J, A < 1). The phe- 
nomenon is essentially driven by d = 1 physics due to the 
extreme sensitivity of a d = 1 system towards a 2fcp per- 
turbation. This is manifest in the superquadratic growth 
of the energy gain AE oc — <5 4 / 3 with the dimerization 
5:=(1-A)/(1+A) §. 

It is the aim of the present work to show that CuGeC>3 
is only approximately a d = 1 compound. On the basis 
of the dispersion perpendicular to the chains this has 
been noted by several authors (e.g. @-§j§J). What has 
passed unnoticed so far is that the spin-spin couplings 
perpendicular to the chains, which are formed along the 
c-axis, influences the dispersion along the chains in a non- 
trivial, non-constant way. 

So far, the deviations of the behavior of CuGe03 
from the physics of simple nearest-neighbor (NN) cou- 
pled Heisenberg chains plus phonons were attributed 
to an intrachain frustration, i.e. an additional next- 
nearest- neighbor (NNN), antiferromagnetic coupling Ji 
@ PlR- But no consensus has been reached which 
value the relative frustration a = J%/J [fi~4| applies 
for CuGe03. The values range from a = 0.36 O] to 
a = 0.17 [[land the values for J from 10.4 meV || to 
21.5 meV |15|. This poses a problem since the NNN cou- 
pling represents a marginal perturbation which engen- 
ders above the critical value a c = 0.241 fill a gap in the 



spin system itself without alternating coupling constants. 
The high temperature data for the magnetic susceptib- 
lity x(T) Jl],[l7j suggests values above a c pi|,p6[ whereas 
there is no indication for a spin gap without coupling 
alternation |l2| ]. An analysis of the dispersion uj(q c ) at 
small q c suggests a = 0.17<a c f§' 

The dispersion in b-direction is roughly 0.3 of the one 
in c-direction; the dispersion in a-direction is roughly 
1/50. First, we focus on the situation in the bc-plane. 
Up to now, the complete dispersions to(q) is interpreted, 
at least qualitatively, as being a sum of two indepen- 
dent contributions ui(q) ss uJbilb) + w c (g c ) leading to the 
schemes in 18 (inset fig. 2) and in Q (fig. 9). 

On the other hand, elastic neutron scattering |l9 20 
provides evidence that the dimerization pattern is shifted 
by one spacing from one chain to the other, i.e. adjacent 
copper ions move in opposite directions [ fl9| |. A segment 
of the resulting dimer distribution is shown in fig. 1. 
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FIG. 1. Segment of the dimerized lattice in the bc-plane. 
The largest coupling is set to unity and the weaker couplings 
are parametrized as indicated. 

The bold lines represent the dimers (J is set to unity); 
the other lines (thin, dashed, dotted) represent the resid- 
ual, weaker couplings (A,/i < 1). The line A is the 
weak bond in the dimerized chain, /i is the interchain 
coupling, fi/3 is a possible NNN hopping, and \a is the 
intrachain frustration discussed before. The microscopic 
superexchange paths for [i and are Cu-02-02-Cu in 
b-direction, see [g(J n S- 1- Since from one of the 02 two 
Cu with different c-coordinate can be reached, a direct 
super-exchange (jx) or a super-exchange with a shift along 
the c-axis (/x/3) is possible |20|. One expects that (3 » 0.5 
because there are two paths for the direct super-exchange 
and one for each shifted super-exchange. 

The analysis starts from the limit of strong dimeriza- 
tion. The experimental systems are closer to S = than 
to S = 1. But the 5=1 case is a good starting point 
due to its simplicity which allows a systematic perturba- 
tive approach, in principle possible down to arbitrarily 
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small S > 0. In practice a third order calculation yields 
reasonable results down to 5 w 0.1 (see fig. 1 in |2l]]). 

At S = 1, the ground state is a product of local sin- 
glets on the dimers and a magnon is a degenerate lo- 
cal single triplet. These triplets acquire a dispersion 
due to the residual couplings which lift the degener- 
acy. To first order one finds a hopping from dimer to 
dimer along the chains of tio = — A(l — 2a) /A |2l[] and of 
tn = — — 2/3) / 4 to the nearest dimers in b-direction. 
To obtain the corresponding dispersion one may view the 
dimers as sites of the effective lattice in fig. 2. 




FIG. 2. Effective lattice if the dimers are taken as elemen- 
tary sites. The hopping elements refer to hops starting at the 
site to- 

Dashed lines correspond to the processes oc fi and 
solid lines to processes oc A. Starting from to onc 
reaches the sites tio, tn and their symmetry equiva- 
lents (reflection at the c-axis and the b-axis) with the 
amplitudes given above. A natural unit cell is spanned 
by e x and e y as indicated. The dashed lines form a 
square lattice to which the solid lines add a compo- 
nent along one of the diagonals. The dispersion reads 
Lu(q) = 2tn(cos(q x ) + cos(q y )) + 2t w cos(q x - q y ) where 
q x / y are the components along the dual basis vectors be- 
longing to e x and to e y , respectively. Changing to the 
basis dual to {e;,, e c } by q x / y = q b /2 ± q c one obtains 
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= 1 + 2* 



ho cos(2<? c ) + Atn cos(q b /2) cos(q c ) . (1) 

Eq. (jl]) illustrates already the main points of this pa- 
per. A dispersion of type (|l|) has minima at q — (0, 0, 0) 
and at (0,1,1/2). Saddle points are at (0,1,0) and at 
(0,0,1/2). These features stay valid in higher pertur- 
bative orders as well. They change the picture used 
so far (cf. fig. 9 in Q) where the symmetry relations 
(0,0,0) ^ (0,0,1/2) and (0,1,0) (0,1,1/2) held. The 
fact that the system is actually 2d (or 3d, see below) does 
not only lower the dispersion along c by some constant 
for constant qb but adds another q c dependence. The 
symmetry about q c — 7r, which one has in d = 1 due 
to dimerization plj , is broken by the term proportional 
to cos(g c ). Since tn is « 1/7 of t w the effect is not 
very large, but visible. Physically this mixture of mo- 
tion in the b- and c-direction is understandable in fig. 1. 



Any hop along b implies a shift also along c because the 
dimers alternate. In addition, the motion along b lowers 
the observable gap. Hence a given gap corresponds to a 
larger dimerization than in d — 1. 

For a more quantitative comparison of theory and ex- 
periment we extend (0) to third order by standard de- 
generate perturbation theory 

t = 1 - (4 - 3a 2 )A 2 /16 - (4 - 3/3 2 )^ 2 /8 - (8 - 8a 

- 6a 2 + 3a 3 )A 3 /64 - (8 - 8/3 - 6/3 2 + 3/3 3 )^ 3 /32 
+ (8a - 16/3 + 18a/3 + 9/3 2 - 12a/3 2 )A^ 2 /32 

+ A 4 /128- 3AV/16- 3/Z/64 (2a) 
tio = -aA/4 - A 2 /8 - /3V/16 - (4 - 4a - a 3 )A 3 /64 

- (2 - /3)/3V/32 + (6a - 6/3 + /3 2 + 2a/3 2 )A^ 2 /32 
+ A 4 /128 + A 2 /i 2 /32 + 3^/128 (2b) 

tn = -/3/i/4 - a/3A/Vl6 - ^ 2 /8 + (4a + 12/3 - 4a/3 
+ 7a 2 /3)A 2 /i/128 + (a - (3)\fJ?/8 - (8 - 16/3 

- 7/3 3 V 3 /128 + 5AV/128 + 3^/128 (2c) 
t 20 = -a 2 X 2 /32 - (2 - a)a 2 A 3 /64 

- 3a/3 2 A^ 2 /64 + 3A 4 /256 (2d) 
t 2 i = -a/3A///16 - (8a -8(3 + 8a(3 + 5a 2 /3)/A 2 ^/256 

+ (2a - 2/3 - 2a(3 + a/3 2 )A/i 2 /64 - 3/3V/256 

+ 3AV/128 (2e) 
t 2 2 = -/3V/32 - 3a/3 2 A^ 2 /64 - (2 - /3)/3V/64 

+ 3/Z/256 (2f) 
*23 = -/3VA6 - 3a/3 2 A^ 2 /64 + (4 - 2(3 + /3 2 )/3/i 3 /32 

+ 3^ 4 /128 (2g) 

t 30 = -a 3 A 3 /128 t 3 i = -3a 2 /3A 2 Ai/128 (2h) 

t 3 2 = -3aP 2 Xfi 2 /128 t 33 = -P 3 fi 3 /128 (2i) 



t 34 = -3/3V/128 



(2j) 



All frustration terms have been included with a — 1 — 2a 
and (3 = 1 — 2(3. This is not too difficult once one has 
the perturbation terms for a = 0, (3 = because the 
frustration terms link the same dimers and they have the 
same effects on them except for the matrix elements. The 
dispersion reads 



J 



■■ [t + 2t w cos(2g c ) + 2t 20 cos(4ij c ) + 2t 30 cos(6g c ) 

-4cos(q 6 /2)(tiicos(g c )+t 2 icos(3g c )+t 3 icos(5<2' c )) 

■ 2 cos(q b )(t 23 + 2t 22 cos(2(7 c ) + 2t 32 cos(4g c )) 

■ 4cos(3<#,/2)(t 34 cos((j c ) + t 33 cos(3g c ))]. (3) 



For jti = 0, a = 0, (^J), (g) can be checked with the 
result for a dimerized NN chain |f22]| . For /i = 0, Aa = 
J/J± , A = one retrieves the result for a NN spin ladder 
|23| . The upper chain is made up by the first sites of 
the dimers; the lower chain of the second sites. For /i = 
0, Aa = J/Jj_ and arbitrary A, a general zigzag chain is 
described. This should give an excellent description of 
the magnon dispersion in Cu 2 (C5Hi 2 N 2 ) 2 Cl4 p4j. 
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Another special case is ((a = 1/2 or A = 0) and (/3 = 
1/2 or li = 0)) where brackets, 'ands', and 'ors' are logical 
signs. The singlet product state is the ground state in the 
d = 1 case \i = 0, A < 1 [E5[. The same can be proven for 
the d — 2 lattice in fig. 2ifa = /3=l/2. For this case the 
fourth order terms are included in (||). It is reasonable 
to include the fourth order terms if a, f3 are close to 1/2. 

20 
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FIG. 3. a: Dispersion along the c-axis resulting from (g) 
and (|) at (J, A, a, fi, /3) = (ll.OmeV, 0.78, 0, 0.30, 0.3) with 
experimental data after fig. 8. Wave vectors in reciprocal 
lattice units (1 r.l.u. = 2n/ lattice spacing), b: Same as in 
fig. 3a along the b-axis (after |J, fig. 5). 

In fig. 3, a comparison between the theory and INS 
results is made. The best fit is obtained for a = al- 
though recent excellent Y(T)-fi ts indicate that a = 0.35 
for T > 30 K is realistic [£6Lp7| . The discrepancy can be 
due to the truncation of the perturbation series at third 
order or to a change of a across the spin-Peierls tran- 
sition induced by the static distortion or induced by a 
renormalization because of non-adiabatic phonon effects. 
This remains to be clarified. 

For the a = fit the terms oc A 2 /i 2 ,A 4 in (g) are 
dropped. The value A = 0.78 corresponding to 5 = 0.12 
is not as small as assumed previously (5 = 0.03 . . .0.05 
PI [l3|,p[) but is of the order of 'first principles' values 
S = 0.21 @ or S ps 0.1 @. Note that the inclusion of 
higher order perturbation terms, enhancing w max — w m in 
at constant A, will lower the appropriate fit value for A 
even further. 



The dispersion along the a-axis is accounted for 
by adding to uj(q) as given by (|J) a term u> a (q) = 
4i Q cos(g a ) cos(g c ) with 4i a = — 0.22mcV M. This term is 
analogous to the third term in (Q) since the dimerization 
pattern alternates also in a-direction [M. 

The agreement is excellent in fig. 3a for q c < 0.3, where 
most data points are available, and it is still very good 
for higher values. The data points show clearly a small 
shift to higher values of q c which results from the terms 
proportional to cos(nq c ) with n odd. Such terms are a 
signature of interchain hopping with alternating dimer- 
ization. The agreement in fig. 3b is also good. 

The fit in fig. 3 corroborates the approach to de- 
scribe the dispersion as hopping on the lattice in fig. 2. 
The parameters, however, (J, A, a, la, /3) are not unam- 
biguously fixed by the dispersion data. Good fits can 
also be obtained for smaller values of li and /3, e.g. 
(J, A, a, n, (3) = (10.7meV, 0.78, 0,0.14,0). This can be 
understood by the fact that the main hopping effect along 
b is cx /Lt(l — 2/3). To fix the value of /3 we deduced the 
coefficients xi ~ —7.8 ± 0.9meV and \2 ~ 26 ± 7meV 2 
in a large T-expansion of the magnetic susceptibility 
X (T) = ( X o/T)(l + Xi/T + X2IT 2 + ...) from the data 
in M. Performing this expansion for the undimerized 
counterpart of the lattice in fig. 1 yields 



Xi = -(Jx + J„ + J a + 2J p )/2 

Xi = {J} + J\(J? + J a ) + 2J (3 (J Q + J M + J A ))/2 



(4a) 
(4b) 



where J\ 



(1 + A) J/2, J a = JXa, J M 



J/i, and 

Jp = Jiifl. For parameters as in fig. 3, one has 
Xi = — 7.5meV, X2 = 30meV 2 while for (3 — one has 
Xi = — 5.4meV, %2 = 6.9mcV 2 . This points to a certain 
amount of interchain frustration (3 not too far away from 
/3 ps 0.5. 

Fig. 4 displays the overall form of uj(q) differing from 
the schemes in JT^Q: a) In fig. 4 to(q) is roughly con- 
stant for q c — 0.25, i.e. close to the maximum value of 
Lo(q), whereas it displays the same dependence on qi, at 
all q c in Jl8|fl. b) According to |l|||] the smaller gap 
E g i PS 2.1mcV is found at q= (0, 1, 1/2) and (0, 1,0); the 
larger gap E g 2 ps 5.5meV is found at q = (0, 0, 1/2) and 
(0,0,0). In fig. 4, E g i and E g2 are swapped at q c = 0. 
This is very important in the interpretation of experi- 
mental data (INS, Raman, ESR, FIR). The magnon DOS 
deduced from fig. 4 is constant at the lower band edge 
(2d minimum), it has a logarithmic singularity at E g 2 (2d 
saddle point), but ressembles to a Id, l/y/uj max — to DOS 
at the maximum due to the approximate degeneracy at 
q c = 0.25. 

Magnetic resonance data [^9 30 shows an excitation 
from the ground state at E g i which was interpreted 
to be the gap at (0,0,0). But a vertical transition at 
q = (0, 0, 0) from 5 = (ground state) to S = 1 (excited 
states) is actually impossible due to total spin conserva- 
tion. 
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the DFG. 



q b /[r.l.u.] q c /[r.l.u. 

FIG. 4. 3D representation of the dispersion u(q) for q a — 
(parameters as in fig. 3). 

For this reason and in view of fig. 4, a re-interpretation 
of this resonance is proposed. In analogy to the magnetic 
resonances due to transitions ground state — > excited 
state in NENP and NINO (see e.g. flpJl) the transi- 
tion reflects the situation at (0, 1, 0) [g3[ because of the 
presence of staggered magnetic fields due to differently 
orientated g-tensors on neighboring chains in b-direction 
||34j| . The Cu-atoms in each unit cell are inequivalent 
since their crystal environments are differently orientated 
{H (cf. the situation for NENP/NINO This ar- 

gument makes the magnetic resonance data compatible 
with fig. 4 (within 10% accuracy as for NENP §j§). 

Further INS investigation are desirable to make a more 
precise analysis possible. It should be noted that the 
dispersion along the a-axis leads to four different gaps 
in fig. 4. For the parameter set used in figs. 3 and 4, 
one has w((0, 1,1/2)) = 2.1meV, w((0,0,0)) = 2.5meV, 
w((0, 0, 1/2)) = 5.8meV, and w((0, 1, 0)) = 6.2meV. 

The magnon dispersion described will, at least qual- 
itatively, also apply to the currently investigated a'- 
NaV2C>5 where the dimerization is equally alternating 
from chain to chain Such an alternation appears to 
be generic for inorganic substances where bond angles 
are more easily changed than bond lengths. 

In summary, we showed that the lattice distortion in 
the spin-Peierls phase causes a dispersion which corre- 
sponds to hopping on a lattice as in fig. 2. Thus the gap 
positions in the Brillouin zone are revisited. A general 
perturbative result is given to enable a quantitative anal- 
ysis of the dispersion. A larger dimerization (8 w 0.1) in 
agreement with orbital calculations is found. This result 
covers also chains and certain ladders which are a very 
active field as well [Q. The magnetic resonance results 
are reinterpreted. 
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